BFF and cellhashR: analysis tools for accurate demultiplexing of cell hashing data.

MOTIVATION: Single-cell sequencing methods provide previously impossible resolution into the transcriptome of individual cells. Cell hashing reduces single-cell sequencing costs by increasing capacity on droplet-based platforms. Cell hashing methods rely on demultiplexing algorithms to accurately classify droplets; however, assumptions underlying these algorithms limit accuracy of demultiplexing, ultimately impacting the quality of single-cell sequencing analyses. RESULTS: We present Bimodal Flexible Fitting (BFF) demultiplexing algorithms BFFcluster and BFFraw, a novel class of algorithms that rely on the single inviolable assumption that barcode count distributions are bimodal. We integrated these and other algorithms into cellhashR, a new R package that provides integrated QC and a single command to execute and compare multiple demultiplexing algorithms. We demonstrate that BFFcluster demultiplexing is both tunable and insensitive to issues with poorly behaved data that can confound other algorithms. Using two well-characterized reference datasets, we demonstrate that demultiplexing with BFF algorithms is accurate and consistent for both well-behaved and poorly behaved input data. AVAILABILITY AND IMPLEMENTATION: cellhashR is available as an R package at https://github.com/BimberLab/cellhashR. cellhashR version 1.0.3 was used for the analyses in this manuscript and is archived on Zenodo at https://www.doi.org/10.5281/zenodo.6402477. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A mechanistic model of PCR for accurate quantification of quantitative PCR data.

BACKGROUND: Quantitative PCR (qPCR) is a workhorse laboratory technique for measuring the concentration of a target DNA sequence with high accuracy over a wide dynamic range. The gold standard method for estimating DNA concentrations via qPCR is quantification cycle () standard curve quantification, which requires the time- and labor-intensive construction of a standard curve. In theory, the shape of a qPCR data curve can be used to directly quantify DNA concentration by fitting a model to data; however, current empirical model-based quantification methods are not as reliable as standard curve quantification. PRINCIPAL FINDINGS: We have developed a two-parameter mass action kinetic model of PCR (MAK2) that can be fitted to qPCR data in order to quantify target concentration from a single qPCR assay. To compare the accuracy of MAK2-fitting to other qPCR quantification methods, we have applied quantification methods to qPCR dilution series data generated in three independent laboratories using different target sequences. Quantification accuracy was assessed by analyzing the reliability of concentration predictions for targets at known concentrations. Our results indicate that quantification by MAK2-fitting is as reliable as standard curve quantification for a variety of DNA targets and a wide range of concentrations. SIGNIFICANCE: We anticipate that MAK2 quantification will have a profound effect on the way qPCR experiments are designed and analyzed. In particular, MAK2 enables accurate quantification of portable qPCR assays with limited sample throughput, where construction of a standard curve is impractical.

Automated identification of axonal growth cones in time-lapse image sequences.

The isolation and purification of axon guidance molecules has enabled in vitro studies of the effects of axon guidance molecule gradients on numerous neuronal cell types. In a typical experiment, cultured neurons are exposed to a chemotactic gradient and their growth is recorded by manual identification of the axon tip position from two or more micrographs. Detailed and statistically valid quantification of axon growth requires evaluation of a large number of neurons at closely spaced time points (e.g. using a time-lapse microscopy setup). However, manual tracing becomes increasingly impractical for recording axon growth as the number of time points and/or neurons increases. We present a software tool that automatically identifies and records the axon tip position in each phase-contrast image of a time-lapse series with minimal user involvement. The software outputs several quantitative measures of axon growth, and allows users to develop custom measurements. For, example analysis of growth velocity for a dissociated E13 mouse cortical neuron revealed frequent extension and retraction events with an average growth velocity of 0.05 +/- 0.14 microm/min. Comparison of software-identified axon tip positions with manually identified axon tip positions shows that the software's performance is indistinguishable from that of skilled human users.